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We measure the dissipation and frequency shift of a magnetically coupled cantilever in the vicinity 
of a silicon chip, down to 25 mK. The dissipation and frequency shift originates from the interaction 
with the unpaired electrons, associated with the dangling bonds in the native oxide layer of the 
silicon, which form a two dimensional system of electron spins. We approach the sample with a 3.43 
/rm-diameter magnetic particle attached to an ultrasoft cantilever, and measure the frequency shift 
and quality factor as a function of temperature and the distance. Using a recent theoretical analysis 
[J. M. de Voogd et ah, arXiv: 1508.07972 (2015)] of the dynamics of a system consisting of a spin and 
a magnetic resonator, we are able to fit the data and extract the relaxation time Ti = 0.39 ± 0.08 
ms and spin density a = 0.14 ±0.01 spins per nm^. Our analysis shows that at temperatures < 500 
mK magnetic dissipation is an important source of non-contact friction. 


Understanding the dissipation and frequency shifts in 
magnetic force experiments is crucial for the develop¬ 
ment of magnetic imaging techniques, e.g. Magnetic 
Resonance Force Microscopy (MRFM). The sensitivity 
of such techniques depends on the friction of the can¬ 
tilevers, which therefore has increased the interest in 
high-quality cantilevers with quality factors exceeding a 
milliorP. However, the quality factor reduces due to non- 
contact friction with the scanned sample which is ex¬ 
plained by dielectric fluctuation^. Far from the surface, 
magnetic dissipation from paramagnetic spins or nano¬ 
magnets on the cantilever ha ve been observed to have a 
large effect on the frictiorPEl. Our report quantitatively 
analyzes the magnetic dissipation of a cantilever in the 
vicinity of a silicon chip, showing that this is the most 
significant non-contact friction at low temperatures for a 
magnet on cantilever geometry. 

Magnetic Force Microscopy (MFM) measures the 
forces resulting from stray fields of a sample that is be¬ 
ing scanned. The coupling of the tip with the magnetic 
field manifests itself as a shift in the resonance frequency 
of the cantilever and as additional dissipation which re¬ 
duces its quality factor Q. For magnetic moments that do 
not change due to the magnetic tip itself, the frequency 
shifts are well understood. However, a more complicated 
model is required when the spins in the sample are para¬ 
magnetic, because the motion of the tip changes the di¬ 
rection of their magnetic moment^. 

In this paper, we show frequency shifts and dissipa¬ 
tion resulting from the dangling electron bonds at the 
surface of a silicon substrate. We are able to extract 
the relaxation time Ti of the electron spins, without us¬ 
ing electron spin resonance techniques. Furthermore, we 
use our analysis to calculate the maximum possible dissi¬ 
pation for a state-of-the-art MRFM setup and diamond 
cantilever. We show that magnetic dissipation can cause 
a drop in quality factor, thereby decreasing the sensitivity 
of an MRFM experiment. We calculate that this dissi¬ 
pation is suppressed when using large external magnetic 
fields at low temperatures. 



FIG. 1. (a) Scanning Electron Microscope image of the mag¬ 

netic particle after it is glued to the cantilever, (b) Optical 
microscope image of the detection chip. The cantilever is po¬ 
sitioned above the center of the pickup coil (•). The pickup 
coil is used for SQUID based detection of the cantilever’s mo¬ 
tion. The vertical wire (dotted arrow) and the copper sample 
(*) are used in other experiments, (c) Sketch of the setup. 
The height is measured from the bottom of the magnetic par¬ 
ticle, which has a diameter of 3.43 pm. (d) The coupling with 
the pickup coil as function of the x-position of the cantilever. 
The red solid line is the calculated flux change in a square 
loop due to a magnetic dipole p on a moving resonator. The 
maximum (scaled to 1) of the curve is at the center of the 
pick-up coil, which can be determined with pm precision. 


In our experiment, a magnet attached to a cantilever 
(Fig. [^) couples via its magnetic field B{r) to magnetic 
moments /x originating from localized electron spins with 
near-negligible interactions. The coupling with a single 
spin can be associated with a stiffness fcg, which results 
in a shift A/ of the natural resonance frequency /o of the 
cantilever, according to A/ = Hj/o, with ko the natural 
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FIG. 2. Frequency sweeps of the cantilever at a temperature 
of 70 mK. When moving towards the sample, the resonance 
frequency fr increases, while Q decreases due to an increasing 
coupling with the surface electron spins. We extract fr and 
Q by fitting the data to a Lorentzian (red solid line). 


stiffness of the cantilever. 

Commonly, the analysis of magnetic interactiorP be¬ 
gins with the interaction energy E = —fi ■ B(r). And 
one calculates the force and stiffness acting on the can¬ 
tilever by taking the first and second derivative with re¬ 
spect to X, the direction of the fundamental mode of the 
cantilever. Assuming that fj, is fixed by a large external 
field, one obtains in this approach a stiffness in the form 
of As = /2 • ■ 

A recent detailed analysis by De Voogd et al.^ which 
starts with the Lagrangian of the full system, taking into 
account the spin’s dynamics as well as the influence of 
the mechanical resonator on the spin, suggests that the 
commonly employed model is not the correct approach 
for paramagnetic spins. For paramagnetic spins, the re¬ 
laxation and the exact dynamics of the spin in the can¬ 
tilever’s magnetic field determine the frequency shifts and 
dissipation. In the case of a two-dimensional system of 
paramagnetic spins, uniformly distributed over an infi¬ 
nite surface, the frequency shift A/ and shift in the in¬ 
verse quality factor A A can be written as: 
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Where T is the temperature, ks is the Boltzmann con¬ 
stant and Ti is the spin’s longitudinal relaxation time. 


The integral is performed over the infinite surface as¬ 
suming a constant spin density cr. We have assumed 
A/ /o, Q ^ 1, and that the inverse of the transverse 
relaxation time is much smaller than the Larmor fre¬ 
quency, which is already the case when is larger than 
1 /xs. 

In this paper, we detect the dangling bonds that are 
present on the surface of a silicon substrate of the detec¬ 
tion chip using MFM down to 25 mK. We use a commer¬ 
cial cryogen-free dilution refrigerator, in which we im¬ 
plemented several vibration isolation measureP. We are 
able to coarse approach towards the sample in three di¬ 
mensions, with a range of 1 mm in x, y and z. For this 
we employ three ‘PiezoKnobs’, from Janssen Precision 
Engineering B.V., while reading out the position using 
three capacitive sensors. 

The cantilever is a silicon micro-machined IBM-type 
with length, width a nd t hickness of 145 /xm, 5 /xm and 
100 nm, respectiveljI^E^. The magnetic particle is a 
spherical particle from a commercial neodymium-alloy 
powdei!^^. We used platinum electron beam induced de¬ 
position using an in-house developed nanomanipulatoili^ 
in a Scanning Electron Microscope (SEM) to attach the 
small magnetic particle on the free end of the cantilever 
and measured the diameter to be 3.43 /xm (Fig. &)■ Sub¬ 
sequently, we magnetized the magnet in the x-direction 
at room temperature in a field of 5 T. 

The readout of the cantilever’s motion is based on a 
Superconducting Quantum Interference Device (SQUID) 
which enables low temperature experimentJ^^. Where in 
conventional MFM setups a laser is used to readout the 
motion, our method is based on the motion of the mag¬ 
netic particle in the vicinity of a small superconducting 
‘pickup’ coil, giving a flux change whenever the cantilever 
moves (Fig. This signal is transformed by an on-chip 
transformer, which matches the pickup coil inductance 
to the high SQUID input inductance. The measnred flux 
noise in the complete setup is less than 4 
where $o is the flux quantum. 

The substrate is high resistivity (> 1 kUcm) (100)- 
oriented n-type (phosphorus doped) silicon. The sub¬ 
strate is cleaned with acetone and DI water, which leaves 
an interface of silicon with its native oxide. To create the 
superconducting structures on the chip, NbTiN is grown 
on the silicon substrate with a thickness of roughly 300 
nm. Patterning is done using standard nano-lithographic 
techniques and reactive ion etching in a SF 6 /O 2 plasma. 
For future MRFM experiments, we added a wire for 
radio-frequency currents and a 300 nm thick copper layer 
capped with gold. The copper is connected via golden 
wire bonds to the sample holder, which itself is connected 
via a silver welded wire to the mixing chamber, ensuring 
good thermalization of the sample. Figure shows an 
optical microscope image of the obtained structure. 

We drive the cantilever using a small piezo element 
glued to the cantilever holder. We sweep the drive fre¬ 
quency using a function generator around the resonance 
frequency fr while measuring the SQUID’s response us- 
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FIG. 3. Resonance frequency fr and quality factor Q versus 
temperature for different heights of the cantilever with respect 
to the sample. For the quality factor, the error bars indicate 
the 95 % confidence intervals of the Lorentzian ht. For the 
frequencies the average error was 0.01 Hz, which is smaller 
than the point size, except for one data point. The solid lines 
are fits to the data with the spin density a, spin relaxation 
time Ti and frequency offset /o as fitting parameters, fr 
and Q are simultaneously fitted for each height. The results 
of the fit can be found in table HI The dashed line is the 
frequency shift calculated with the commonly used expression 
ks = pfi ■ , with p = tanh ^ ^ and with a ten times 

smaller than we find in our analysis. 


ing a Lock-In amplifier. We fitted the square of the 
squid’s signal with a Lorentzian curve in order to ex¬ 
tract fr and Q. The amplitude of the Lorentzian is de¬ 
termined by the coupling between the magnet and the 
pickup coil, which is proportional to the energy coupling, 
and can be used to determine the position of the can¬ 
tilever by scanning the cantilever in the xy-plane, see 
figure [^. 

For the experiment presented in this paper, we posi¬ 
tioned the cantilever above the center of the pickup coil, 
to minimize possible repulsive forces from the supercon¬ 
ducting wires. By gently decreasing the height of the 
cantilever until the signal is completely lost, we deter¬ 
mine the relative height of the magnetic particle with 
respect to the surface. The sample holder is placed on 
a finestage, machined out of aluminum, which can be 
moved in all spatial directions by actuating laminated 


piezoelectric extension stacks. Using this, we can now 
have good control of the height up to the full range of 
the finestage of 2.3 /rirfSI. 


We swept the drive frequency at a drive amplitude 
small enough to avoid non-linear responses of the can¬ 
tilever’s motion, while measuring the SQUID signal. We 
measured with a sampling time of 2 s every 0.02 Hz. 
Fitting the data with a Lorentzian, we obtain fr and 
Q = ■ At each height, the temperature was var¬ 

ied from the lowest achievable temperature « 25 mK, up 
to 500 mK. Above 500 mK, the aluminum shielding of 
the experiment starts to become non-superconducting. 
An example of the data with the Lorentzian hts at all 
used heights at 70 mK is shown in figure 

The results of our measurements described above are 
shown in figure [3| together with the fits according to equa¬ 
tions Q and (0. At every height z and temperature T 
we calculate the value for C according to equation 
The quality factor far from the surface Qq = 2.8 • 1^. 


The stiffness kg = 7.0 • 10“® Nm“^ of the cantilever is 
calculated using ko = nieff (27r/o)^ with /o = 3.0 kHz 
and rueff = 2.0 • 10“^^ kg. The effective mass We// is 
calculated using the geometry of the cantilever and the 
magnetic particle. The magnetic particle is taken as a 
spherical dipole with magnetic moment m. According 
to the model, the temperature at which the resonance 
frequency close to the sample has a maximum, is inde¬ 
pendent of a and Ti, but is dependent on the absolute 
value of m and the distance to the sample. We find 
m = 1.9 • 10“^^ JT“^. From this we hnd an effective sat¬ 
uration magnetization of 1.15 T for a sphere that is fully 
magnetic. Alternatively we can assume pL^Mgat = 1-3 T 
and an outer layer of 200 nm which is magnetically dead. 
The magnetic moment of the dangling bondJ^^ is equal 
to the Bohr magneton /i = 9.274 • 10“^'^JT“^. 

The solid lines in hgure ^are hts to the data accord¬ 
ing to equations Q and with a, Ti and /o as the 
only htting parameters. All htting parameters are sepa¬ 
rately htted for each height, for both the frequency data 
and the quality factor data, /o is a temperature inde¬ 
pendent parameter different for each height, which we 
attribute to an unknown mechanism, since the coupling 


TABLE I. Obtained values for the spin density cr and relax¬ 
ation time Ti for every height ^ above the sample. See for 
the individual fits hgure. The bottom row shows the average 
value and the standard deviation. 


Height (pm) 

spin density (nm’ 

Relaxation time (ms) 

0.08 

0.142 

0.42 

0.19 

0.137 

0.52 

0.38 

0.140 

0.48 

0.57 

0.142 

0.42 

0.77 

0.136 

0.38 

1.15 

0.130 

0.32 

1.72 

0.133 

0.28 

2.30 

0.168 

0.33 

mean: 

0.14 ± 0.01 

0.39 ±0.08 
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to the SQUID is too small of an effect at these distances 
and has a height dependence with opposite sign to the 
one observed. The results of the fits for Ti and a can be 
found in table |Tj We left a as fitting parameter for each 
height, to verify the correctness of our analysis, since 
this number should be the same for each height. We 
see that Ti slightly increases when the magnetic parti¬ 
cle approached the surface, as is also observed for bulk 
spins in electron spin resonance experimental. Ti could 
depend on temperature, but by taking the ratio of equa¬ 
tion 0 with equation 0 we extract Ti for each measure¬ 
ment, and we find that Ti is constant with temperature 
to within 20%. The average values of all individual fits 
are a = 0.14 ± 0.01 spins per nm^ and Ti = 0.39 ± 0.08 
ms. The found value for a is similar to value s measured 
using Electron Paramagnetic Resonanc d^^hd 

The dashed line in figure]^ is the frequency shift calcu¬ 
lated with the commonly used expression = pfJ, ■ , 

with p = tanh ^ ) • Important is that for this curve, 

the spin density is ten times smaller than we find with 
our analysis. 

The deviation of the data from the fit for low tem¬ 
peratures and small values for z can be understood by 
considering that we do not have only spins at the sur¬ 
face. Electron spins inside the bulk will cause deviations 
to the fits, already when the density is in the order of 
10 '^ spins per fim^ which is less than 1 ppm of the silicon 
atoms. Considering the nuclear spins, the 4.7% natural 
abundance of the ^®Si isotope can only account for less 
than 1 percent deviation. 

Note that in electron spin resonance studies with our 
MRFM setup, a value for Ti in the order of seconds was 
reportecP^. With our new analysis we believe that it 
is possible that the reported long lived frequency shifts 
could be caused by nuclei polarized by interactions with 
these electron spins, and that these electron spins were 
actually much shorter lived, as is reported for nitroxide- 
doped perdeuterated polystyrene filmJi^. 

Our analysis suggests that the spin mediated dissipa¬ 
tion is the main mechanism leading to a significant reduc¬ 
tion in the quality factor of the cantilever. Previous work 
at higher temperature^ reports dielectric fluctuations as 
the main non-contact dissipation mechanism. We do not 
see any evidence in our measurements for this. Possibly, 
the use of a laser in the setup to read out the cantilever 
causes extra charge fluctuations. Furthermore, we work 
at lower temperatures, where the large spin polarization 
enhances the magnetic dissipation and possibly reduce 
fluctuating charges. 

We calculated the magnetic dissipation for a magnetic 
imaging experiments at higher temperature and a dif¬ 
ferent tip-sample geometry. The results can be found 
in table We used the experimental parameters for a 
state-of-the-art MRFIVP^. In this apparatus, the bare 
non-magnetic cantilever is centered approximately 50 nm 
above a magnetic particle on the substrate, which is as¬ 
sumed for simplicity to be a spherical particle with a 


TABLE II. Calculated quality factor Q for three different tem¬ 
peratures and external magnetic fields assuming magnetic dis¬ 
sipation as the only source for non-contact friction. Calcula¬ 
tions are based on a state-of-the art MRFM apparatus with 
a ‘sample on cantilever’ geometrjl^ and a cantilevei^ with an 
internal quality factor Qo = 1.5 • 10®. 


Calculated quality factors (TO®) 




T = 10 mK' 

T = 300 mK 

T = 4.2 K 

Bext 

= t) T 

0.49 

0".20 

0.98 

^ext 

= 0.1 T 

1.50 

0.19 

0.91 

Bext 

= 6 T 

1.50 

1.50 

1.17 


radius of 100 nm. This setup is equivalent to a mag¬ 
netic dipole attached to the cantilever itself approaching 
a surface with the shape of the cantilever. The magnetic 
dipole and external field are oriented in the z-direction 
while the fundamental mode of the cantilever is in the x- 
direction. For the cantilever, we used the parameters of 
a recently developed diamond cantilevei^ which is shown 
to have at low temperatures an intrinsic quality factor 
Qo = 1.5 • 106 , resonance frequency /o = 32 kHz and 
stiffness ko = 6.7 • 10“^ Nm“^. A spin density a — 0.14 
nm“^ is used, which is found in this report to be the 
density for the silicon surface, but it is also close to the 
typical values found for diamond surface^^. Only spins 
at the very end of the cantilever are considered since this 
surface contributes most to the dissipation, which is 0.66 
/im thick and 12 pm wide. Although equation Q cannot 
be used since we do not have a uniform infinite surface 
anymore, according to the original expression^ one can 
continue to use equation ([^ for the dissipation replacing 
the integral in equation (pr over the end of the cantilever. 
The relaxation time is chosen such that the dissipation 
is maximum: Ti = (27r/o) ^ = 5.0 ps. 

The values in table |ll] show that the magnetic dissi¬ 
pation could be an important source of non-contact fric¬ 
tion. Furthermore we see that applying external fields 
can reduce the magnetic dissipation. Considering these 
calculations we believe that the magnetic dissipation we 
find at low temperatures can be avoided with the cor¬ 
rect choice for the substrate and the use of large external 
magnetic fields. 

To summarize, we have shown how the dissipation and 
frequency shift mediated by spins in magnetic force ex¬ 
periments can be fully understood. The new analysis 
suggest that in order to achieve higher sensitivity in mag¬ 
netic imaging techniques, one should not only focus on 
improving the intrinsic losses of the micro-mechanical 
cantilever, but also on the reduction of electron spins in 
the sample. Furthermore we have shown how the spin’s 
relaxation time can be extracted without the use of res¬ 
onance techniques. For silicon substrates with native ox¬ 
ides, we find a relaxation time of Ti — 0.39±0.08 ms and 
a spin density of ct = 0.14 ± 0.01 per nm^. The under¬ 
standing of the spin mediated dissipation is important to 
further improve the mechanical resonators in magnetic 
imaging experiments. 
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